The Evolution, Expression Patterns, and Domestication Selection Analysis of the Annexin Gene Family in the Barley Pan-Genome.

Plant annexins constitute a conserved protein family that plays crucial roles in regulating plant growth and development, as well as in responses to both biotic and abiotic stresses. In this study, a total of 144 annexin genes were identified in the barley pan-genome, comprising 12 reference genomes, including cultivated barley, landraces, and wild barley. Their chromosomal locations, physical-chemical characteristics, gene structures, conserved domains, and subcellular localizations were systematically analyzed to reveal the certain differences between wild and cultivated populations. Through a cis-acting element analysis, co-expression network, and large-scale transcriptome analysis, their involvement in growth, development, and responses to various stressors was highlighted. It is worth noting that HvMOREXann5 is only expressed in pistils and anthers, indicating its crucial role in reproductive development. Based on the resequencing data from 282 barley accessions worldwide, genetic variations in thefamily were investigated, and the results showed that 5 out of the 12 identified HvMOREXanns were affected by selection pressure. Genetic diversity and haplotype frequency showed notable reductions between wild and domesticated barley, suggesting that a genetic bottleneck occurred on the annexin family during the barley domestication process. Finally, qRT-PCR analysis confirmed the up-regulation of HvMOREXann7 under drought stress, along with significant differences between wild accessions and varieties. This study provides some insights into the genome organization and genetic characteristics of the annexin gene family in barley at the pan-genome level, which will contribute to better understanding its evolution and function in barley and other crops.

Genome-wide association studies of peduncle length in wheat under rain-fed and irrigating field conditions.

Drought is one of the most destructive environmental factors limiting wheat production and food security globally. Peduncle length (PLE) is an important morphological trait to determine plant architecture, photosynthate transport, and yield formation, which is also considered a useful index for drought tolerance in wheat. However, the genetic basis of wheat PLE is not well studied at present. Here, a large-scale genome-wide association study (GWAS) of PLE was performed using a panel of 282 wheat accessions with the Wheat 660K SNP array genotyping under rain-fed and irrigating field conditions. Totally, 1,301 significant marker-trait associations (MTAs) were identified using the threshold of p-value < 4.16 × 10-4, five of which were high-confidence. Furthermore, combining GWAS intervals, previously reported QTLs, expression levels, homologous genes, and selected sweep analysis, a total of 5 candidate genes were detected to associate with drought stress. Moreover, the expression levels of TraesCS2A02G082100 were significantly up-regulated under drought conditions and co-localized in the selected sweep region, suggesting it is a drought-responsive gene. Our results shed light on the genetic basis underlying wheat drought tolerance, which accelerates the marker-assistant selection and genetic improvement through genomic breeding in wheat.

Integrate Small RNA and Degradome Sequencing to Reveal Drought Memory Response in Wheat (Triticum aestivum L.).

Drought has gradually become one of the most severe abiotic stresses on plants. Plants that experience stress training can exhibit enhanced stress tolerance. According to MicroRNA (miRNA) sequencing data, this study identified 195 candidate drought memory-related miRNAs in wheat, and targets of 64 (32.8%) candidate miRNAs were validated by degradome sequencing. Several drought memory-related miRNAs such as tae-miR9676-5p, tae-MIR9676-p3_1ss21GA, tae-miR171a, tae-miR531_L-2, tae-miR408_L-1, PC-3p-5049_3565, tae-miR396c-5p, tae-miR9778, tae-miR164a-5p, and tae-miR9662a-3p were validated as having a strong response to drought memory by regulating the expression of their target genes. In addition, overexpression of drought memory-related miRNA, tae-miR531_L-2, can remarkably improve the drought tolerance of transgenic Arabidopsisthaliana. Drought memory can regulate plant cellular signal transduction, plant biosynthetic processes, and other biological processes to cope with drought via transcriptional memory. In addition, drought memory-related miRNAs can promote starch and sucrose catabolism and soluble sugar accumulation and regulate proline homeostasis to improve plant drought resistance. Our results could contribute to an understanding of drought memory in wheat seedlings and may provide a new strategy for drought-resistant breeding.

The Landscape of Autophagy-Related (ATG) Genes and Functional Characterization of TaVAMP727 to Autophagy in Wheat.

Autophagy is an indispensable biological process and plays crucial roles in plant growth and plant responses to both biotic and abiotic stresses. This study systematically identified autophagy-related proteins (ATGs) in wheat and its diploid and tetraploid progenitors and investigated their genomic organization, structure characteristics, expression patterns, genetic variation, and regulation network. We identified a total of 77, 51, 29, and 30 ATGs in wheat, wild emmer, T. urartu and A. tauschii, respectively, and grouped them into 19 subfamilies. We found that these autophagy-related genes (ATGs) suffered various degrees of selection during the wheat's domestication and breeding processes. The genetic variations in the promoter region of Ta2A_ATG8a were associated with differences in seed size, which might be artificially selected for during the domestication process of tetraploid wheat. Overexpression of TaVAMP727 improved the cold, drought, and salt stresses resistance of the transgenic Arabidopsis and wheat. It also promoted wheat heading by regulating the expression of most ATGs. Our findings demonstrate how ATGs regulate wheat plant development and improve abiotic stress resistance. The results presented here provide the basis for wheat breeding programs for selecting varieties of higher yield which are capable of growing in colder, drier, and saltier areas.

A large-scale genomic association analysis identifies the candidate causal genes conferring stripe rust resistance under multiple field environments.

The incorporation of resistance genes into wheat commercial varieties is the ideal strategy to combat stripe or yellow rust (YR). In a search for novel resistance genes, we performed a large-scale genomic association analysis with high-density 660K single nucleotide polymorphism (SNP) arrays to determine the genetic components of YR resistance in 411 spring wheat lines. Following quality control, 371 972 SNPs were screened, covering over 50% of the high-confidence annotated gene space. Nineteen stable genomic regions harbouring 292 significant SNPs were associated with adult-plant YR resistance across nine environments. Of these, 14 SNPs were localized in the proximity of known loci widely used in breeding. Obvious candidate SNP variants were identified in certain confidence intervals, such as the cloned gene Yr18 and the major locus on chromosome 2BL, despite a large extent of linkage disequilibrium. The number of causal SNP variants was refined using an independent validation panel and consideration of the estimated functional importance of each nucleotide polymorphism. Interestingly, four natural polymorphisms causing amino acid changes in the gene TraesCS2B01G513100 that encodes a serine/threonine protein kinase (STPK) were significantly involved in YR responses. Gene expression and mutation analysis confirmed that STPK played an important role in YR resistance. PCR markers were developed to identify the favourable TraesCS2B01G513100 haplotype for marker-assisted breeding. These results demonstrate that high-resolution SNP-based GWAS enables the rapid identification of putative resistance genes and can be used to improve the efficiency of marker-assisted selection in wheat disease resistance breeding.

Genomic and Phenotypic Divergence in Wild Barley Driven by Microgeographic Adaptation.

Microgeographic adaptation is a fundamental driving force of evolution, but the underlying causes remain undetermined. Here, the phenotypic, genomic and transcriptomic variations of two wild barley populations collected from sharply divergent and adjacent micro-geographic sites to identify candidate genes associated with edaphic local adaptation are investigated. Common garden and reciprocal transplant studies show that large phenotypic differentiation and local adaptation to soils occur between these populations. Genetic, phylogenetic and admixture analyses based on population resequencing show that significant genetic divergences occur between basalt and chalk populations. These divergences are consistent with the phenotypic variations observed in the field. Genome sweep analyses reveal 162.7 Mb of selected regions driven by edaphic local adaptation, in which 445 genes identified, including genes associated with root architecture, metal transport/detoxification, and ABA signaling. When the phenotypic, genomic and transcriptomic data are combined, HvMOR, encoding an LBD transcription factor, is determined to be the vital candidate for regulating the root architecture to adapt to edaphic conditions at the microgeographic scale. This study provides new insights into the genetic basis of edaphic adaptation and demonstrates that edaphic factors may contribute to the evolution and speciation of barley.

The improved assembly of 7DL chromosome provides insight into the structure and evolution of bread wheat.

Wheat is one of the most important staple crops worldwide and also an excellent model species for crop evolution and polyploidization studies. The breakthrough of sequencing the bread wheat genome and progenitor genomes lays the foundation to decipher the complexity of wheat origin and evolutionary process as well as the genetic consequences of polyploidization. In this study, we sequenced 3286 BACs from chromosome 7DL of bread wheat cv. Chinese Spring and integrated the unmapped contigs from IWGSC v1 and available PacBio sequences to close gaps present in the 7DL assembly. In total, 8043 out of 12 825 gaps, representing 3 491 264 bp, were closed. We then used the improved assembly of 7DL to perform comparative genomic analysis of bread wheat (Ta7DL) and its D donor, Aegilops tauschii (At7DL), to identify domestication signatures. Results showed a strong syntenic relationship between Ta7DL and At7DL, although some small rearrangements were detected at the distal regions. A total of 53 genes appear to be lost genes during wheat polyploidization, with 23% (12 genes) as RGA (disease resistance gene analogue). Furthermore, 86 positively selected genes (PSGs) were identified, considered to be domestication-related candidates. Finally, overlapping of QTLs obtained from GWAS analysis and PSGs indicated that TraesCS7D02G321000 may be one of the domestication genes involved in grain morphology. This study provides comparative information on the sequence, structure and organization between bread wheat and Ae. tauschii from the perspective of the 7DL chromosome, which contribute to better understanding of the evolution of wheat, and supports wheat crop improvement.

N6-methyladenosine regulatory machinery in plants: composition, function and evolution.

N6-methyladenosine (m6A) RNA methylation, one of the most pivotal internal modifications of RNA, is a conserved post-transcriptional mechanism to enrich and regulate genetic information in eukaryotes. The scope and function of this modification in plants has been an intense focus of study, especially in model plant systems. The characterization of plant m6A writers, erasers and readers, as well as the elucidation of their functions, is currently one of the most fascinating hotspots in plant biology research. The functional analysis of m6A in plants will be booming in the foreseeable future, which could contribute to crop genetic improvement through epitranscriptome manipulation. In this review, we systematically analysed and summarized recent advances in the understanding of the structure and composition of plant m6A regulatory machinery, and the biological functions of m6A in plant growth, development and stress response. Finally, our analysis showed that the evolutionary relationships between m6A modification components were highly conserved across the plant kingdom.

Conservation analysis of long non-coding RNAs in plants.

Long non-coding RNAs (lncRNAs) are gene regulators that have vital roles in development and adaptation to the environment in eukaryotes. However, the structural and evolutionary analyses of plant lncRNAs are limited. In this study, we performed an analysis of lncRNAs in five monocot and five dicot species. Our results showed that plant lncRNA genes were generally shorter and had fewer exons than protein-coding genes. The numbers of lncRNAs were positively correlated with the numbers of protein-coding genes in different plant species, despite a high range of variation. Sequence conservation analysis showed that the majority of lncRNAs had high sequence conservation at the intra-species and sub-species levels, reminiscent of protein-coding genes. At the inter-species level, a subset of lncRNAs were highly diverged at the nucleotide level, but conserved by position. Interestingly, we found that plant lncRNAs have identical splicing signals, and those which can form precursors or targets of miRNAs have a conservative identity in different species. We also revealed that most of the lowly expressed lncRNAs were tissue-specific, while those highly conserved were constitutively transcribed. Meanwhile, we characterized a subset of rice lncRNAs that were co-expressed with their adjacent protein-coding genes, suggesting they may play cis-regulatory roles. These results will contribute to understanding the biological significance and evolution of lncRNAs in plants.

Characterization of microRNAs and their targets in wild barley (Hordeum vulgare subsp. spontaneum) using deep sequencing.

MicroRNAs (miRNA) are a class of small, endogenous RNAs that play a negative regulatory role in various developmental and metabolic processes of plants. Wild barley (Hordeum vulgare subsp. spontaneum), as the progenitor of cultivated barley (Hordeum vulgare subsp. vulgare), has served as a valuable germplasm resource for barley genetic improvement. To survey miRNAs in wild barley, we sequenced the small RNA library prepared from wild barley using the Illumina deep sequencing technology. A total of 70 known miRNAs and 18 putative novel miRNAs were identified. Sequence analysis revealed that all of the miRNAs identified in wild barley contained the highly conserved hairpin sequences found in barley cultivars. MiRNA target predictions showed that 12 out of 52 miRNA families were predicted to target transcription factors, including 8 highly conserved miRNA families in plants and 4 wheat-barley conserved miRNA families. In addition to transcription factors, other predicted target genes were involved in diverse physiological and metabolic processes and stress defense. Our study for the first time reported the large-scale investigation of small RNAs in wild barley, which will provide essential information for understanding the regulatory role of miRNAs in wild barley and also shed light on future practical utilization of miRNAs for barley improvement.

Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L.

RNA editing is an important way to convert cytidine (C) to uridine (U) at specific sites within RNA molecules at a post-transcriptional level in the chloroplasts of higher plants. Although it has been systematically studied in many plants, little is known about RNA editing in the wheat D genome donor Aegilops tauschii L. Here, we investigated the chloroplast RNA editing of Ae. tauschii and compared it with other wheat relatives to trace the evolution of wheat. Through bioinformatics prediction, a total of 34 C-to-U editing sites were identified, 17 of which were validated using RT-PCR product sequencing. Furthermore, 60 sites were found by the RNA-Seq read mapping approach, 24 of which agreed with the prediction and six were validated experimentally. The editing sites were biased toward tCn or nCa trinucleotides and 5'-pyrimidines, which were consistent with the flanking bases of editing sites of other seed plants. Furthermore, the editing events could result in the alteration of the secondary structures and topologies of the corresponding proteins, suggesting that RNA editing might impact the function of target genes. Finally, comparative analysis found some evolutionarily conserved editing sites in wheat and two species-specific sites were also obtained. This study is the first to report on RNA editing in Aegilops tauschii L, which not only sheds light on the evolution of wheat from the point of view of RNA editing, but also lays a foundation for further studies to identify the mechanisms of C-to-U alterations.

Global Identification of MicroRNAs and Their Targets in Barley under Salinity Stress.

Salinity is a major limiting factor for agricultural production worldwide. A better understanding of the mechanisms of salinity stress response will aid efforts to improve plant salt tolerance. In this study, a combination of small RNA and mRNA degradome sequencing was used to identify salinity responsive-miRNAs and their targets in barley. A total of 152 miRNAs belonging to 126 families were identified, of which 44 were found to be salinity responsive with 30 up-regulated and 25 down-regulated respectively. The majority of the salinity-responsive miRNAs were up-regulated at the 8h time point, while down-regulated at the 3h and 27h time points. The targets of these miRNAs were further detected by degradome sequencing coupled with bioinformatics prediction. Finally, qRT-PCR was used to validate the identified miRNA and their targets. Our study systematically investigated the expression profile of miRNA and their targets in barley during salinity stress phase, which can contribute to understanding how miRNAs respond to salinity stress in barley and other cereal crops.

Preparation of high molecular weight gDNA and bacterial artificial chromosome (BAC) libraries in plants.

Bacterial artificial chromosome (BAC) libraries are extremely valuable large-insert DNA libraries for physical mapping, positional cloning, comparative genomic analysis, complete genome sequencing, and evolutionary studies. Due to their stability and relative simplicity BAC libraries are most preferred over other approaches for cloning large genomic DNA fragments for large-insert libraries. Isolation of intact high molecular weight (HMW) DNA is a critical step underlying the success of large-insert genomic DNA library construction. It requires the isolation of purified nuclei, embedding them into LMP agarose plugs, restriction digestion of the plugs, and quite often size selection using PFGE and electro-elution of insert DNA. The construction of BAC libraries is complex and challenging for most molecular laboratories. To facilitate the construction of BAC libraries, we present a step-by-step protocol for isolation of HMW DNA and construction of plant BAC libraries.

A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor.

The current limitations in genome sequencing technology require the construction of physical maps for high-quality draft sequences of large plant genomes, such as that of Aegilops tauschii, the wheat D-genome progenitor. To construct a physical map of the Ae. tauschii genome, we fingerprinted 461,706 bacterial artificial chromosome clones, assembled contigs, designed a 10K Ae. tauschii Infinium SNP array, constructed a 7,185-marker genetic map, and anchored on the map contigs totaling 4.03 Gb. Using whole genome shotgun reads, we extended the SNP marker sequences and found 17,093 genes and gene fragments. We showed that collinearity of the Ae. tauschii genes with Brachypodium distachyon, rice, and sorghum decreased with phylogenetic distance and that structural genome evolution rates have been high across all investigated lineages in subfamily Pooideae, including that of Brachypodieae. We obtained additional information about the evolution of the seven Triticeae chromosomes from 12 ancestral chromosomes and uncovered a pattern of centromere inactivation accompanying nested chromosome insertions in grasses. We showed that the density of noncollinear genes along the Ae. tauschii chromosomes positively correlates with recombination rates, suggested a cause, and showed that new genes, exemplified by disease resistance genes, are preferentially located in high-recombination chromosome regions.

Dispersion and domestication shaped the genome of bread wheat.

Despite the international significance of wheat, its large and complex genome hinders genome sequencing efforts. To assess the impact of selection on this genome, we have assembled genomic regions representing genes for chromosomes 7A, 7B and 7D. We demonstrate that the dispersion of wheat to new environments has shaped the modern wheat genome. Most genes are conserved between the three homoeologous chromosomes. We found differential gene loss that supports current theories on the evolution of wheat, with greater loss observed in the A and B genomes compared with the D. Analysis of intervarietal polymorphisms identified fewer polymorphisms in the D genome, supporting the hypothesis of early gene flow between the tetraploid and hexaploid. The enrichment for genes on the D genome that confer environmental adaptation may be associated with dispersion following wheat domestication. Our results demonstrate the value of applying next-generation sequencing technologies to assemble gene-rich regions of complex genomes and investigate polyploid genome evolution. We anticipate the genome-wide application of this reduced-complexity syntenic assembly approach will accelerate crop improvement efforts not only in wheat, but also in other polyploid crops of significance.

Development of chromosome-arm-specific microsatellite markers in Triticum aestivum (Poaceae) using NGS technology.

PREMISE OF THE STUDY: The aim of this study was to assess the feasibility of developing chromosome-arm-specific microsatellite markers in wheat on a large scale based on chromosome survey sequences obtained with next-generation sequencing (NGS) technology. METHODS AND RESULTS: The Illumina Hi Seq2000 sequencing platform was used to sequence DNA of isolated wheat chromosome-arm 7DL. The data were assembled and microsatellite loci were identified computationally. In total, 16315 microsatellites were identified from 161061 assembled contigs. Thirty-three markers were randomly selected for validation across 20 diverse wheat cultivars. Two nulli-tetrasomic stocks were also screened to validate the specificity of the newly developed markers. CONCLUSIONS: This is the first study on identification of chromosome-arm-specific microsatellite markers using NGS technology. These new chromosome-arm-specific markers will facilitate saturation of the 7DL genetic map, and their availability will support genetic mapping and positional cloning in wheat.

Complete chloroplast genome sequence of a major invasive species, crofton weed (Ageratina adenophora).

BACKGROUND: Crofton weed (Ageratina adenophora) is one of the most hazardous invasive plant species, which causes serious economic losses and environmental damages worldwide. However, the sequence resource and genome information of A. adenophora are rather limited, making phylogenetic identification and evolutionary studies very difficult. Here, we report the complete sequence of the A. adenophora chloroplast (cp) genome based on Illumina sequencing. METHODOLOGY/PRINCIPAL FINDINGS: The A. adenophora cp genome is 150, 689 bp in length including a small single-copy (SSC) region of 18, 358 bp and a large single-copy (LSC) region of 84, 815 bp separated by a pair of inverted repeats (IRs) of 23, 755 bp. The genome contains 130 unique genes and 18 duplicated in the IR regions, with the gene content and organization similar to other Asteraceae cp genomes. Comparative analysis identified five DNA regions (ndhD-ccsA, psbI-trnS, ndhF-ycf1, ndhI-ndhG and atpA-trnR) containing parsimony-informative characters higher than 2%, which may be potential informative markers for barcoding and phylogenetic analysis. Repeat structure, codon usage and contraction of the IR were also investigated to reveal the pattern of evolution. Phylogenetic analysis demonstrated a sister relationship between A. adenophora and Guizotia abyssinica and supported a monophyly of the Asterales. CONCLUSION: We have assembled and analyzed the chloroplast genome of A. adenophora in this study, which was the first sequenced plastome in the Eupatorieae tribe. The complete chloroplast genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.

Identification and characterization of microRNAs from barley (Hordeum vulgare L.) by high-throughput sequencing.

MicroRNAs (miRNAs) are a class of endogenous RNAs that regulates the gene expression involved in various biological and metabolic processes. Barley is one of the most important cereal crops worldwide and is a model organism for genetic and genomic studies in Triticeae species. However, the miRNA research in barley has lagged behind other model species in grass family. To obtain more information of miRNA genes in barley, we sequenced a small RNA library created from a pool of equal amounts of RNA from four different tissues using Solexa sequencing. In addition to 126 conserved miRNAs (58 families), 133 novel miRNAs belonging to 50 families were identified from this sequence data set. The miRNA* sequences of 15 novel miRNAs were also discovered, suggesting the additional evidence for existence of these miRNAs. qRT-PCR was used to examine the expression pattern of six randomly selected miRNAs. Some miRNAs involved in drought and salt stress response were also identified. Furthermore, the potential targets of these putative miRNAs were predicted using the psRNATarget tools. Our results significantly increased the number of novel miRNAs in barley, which should be useful for further investigation into the biological functions and evolution of miRNAs in barley and other species.

A MAPK gene from Dead Sea fungus confers stress tolerance to lithium salt and freezing-thawing: Prospects for saline agriculture.

The Dead Sea is one of the most saline lakes on earth ( approximately 340 g/liter salinity) and is approximately 10 times saltier than the oceans. Eurotium herbariorum, a common fungal species, was isolated from its water. EhHOG gene, encoding a mitogen-activated protein kinase (MAPK) that plays an essential role in the osmoregulatory pathway in yeast and many other eukaryotes, was isolated from E. herbariorum. The deduced amino acid sequences of EhHOG indicated high similarity with homologous genes from Aspergillus nidulans, Saccharomyces cerevisiae, and Schizosaccharomyces pombe and contained a TGY motif for phosphorylation by MAPK kinase. When EhHOG was expressed in S. cerevisiae hog1Delta mutant, the growth and aberrant morphology of hog1Delta mutant was restored under high osmotic stress condition. Moreover, intracellular glycerol content in the transformant increased to a much higher level than that in the mutant during salt-stress situations. hog1Delta mutant complemented by EhHOG outperformed the wild type or had higher genetic fitness under high Li(+) and freezing-thawing conditions. The present study revealed the putative presence of a high-osmolarity glycerol response (HOG) pathway in E. herbariorum and the significance of EhHOG in osmotic regulation, heat stress, freeze stress, and oxidative stress. The Dead Sea is becoming increasingly more saline while the fungi living in it evolutionarily adapt to its high-saline environment, particularly with the extraordinarily high Li(+) concentration. The Dead Sea is potentially an excellent model for studies of evolution under extreme environments and is an important gene pool for future agricultural genetic engineering prospects.

Wild barley eibi1 mutation identifies a gene essential for leaf water conservation.

Drought is a major abiotic stress that limits plant growth and crop productivity. A spontaneous wilty mutant (eibi1) hypersensitive to drought was identified from wild barley (Hordeum spontaneum Koch). eibi1 showed the highest relative water loss rate among the known wilty mutants, which indicates that eibi1 is one of the most drought-sensitive mutants. eibi1 had the same abscisic acid (ABA) level, the same ability to accumulate stress-induced ABA, and the same stomatal movement in response to light, dark, drought, and exogenous ABA as the wild type, revealing that eibi1 was neither an ABA-deficient nor an ABA-insensitive mutant. The eibi1 leaves had a larger chlorophyll efflux rate in 80% ethanol than the wild-type leaves; and the transpiration rate of eibi1 was more closely related to chlorophyll efflux rate than to stomatal density, demonstrating that the cuticle of eibi1 was defective. eibi1 will be a promising candidate to study the actual barrier layer in the cuticle that limits water loss of the plant. Exogenous ABA reduced leaf length growth in eibi1 more than in the wild type, implying an interaction on plant growth of ABA signal transduction and the eibi1 product. One may infer that the eibi1 product may reverse the growth inhibition induced by ABA.